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Abstract. We present an investigation of velocity fields in early to late M-type hydrodynamic stellar 
atmosphere models. These velocities will be expressed in classical terms of micro- and macro- 
turbulent velocities for usage in ID spectral synthesis. The M-star model parameters range between 
log^ of 3.0 - 5.0 and T^-ff of 2500 K - 4000 K. We characterize the T^g- and logg-dependence of the 
hydrodynamical velocity fields in these models with a binning method, and for the determination 
of micro-turbulent velocities, the Curve of Growth method is used. The macro-turbulent velocities 
are obtained by convolutions with Gaussian profiles. Velocity fields in M-stars strongly depend on 
logg and Teff. Their velocity amplitudes increase with decreasing logg and increasing Teff. The 3D 
hydrodynamical and ID macro-turbulent velocities range from ~ 100 m/s for cool high gravity 
models to ~ 800 m/s-1000 m/s for hot models or models with low logg. The micro-turbulent 
velocities range in the order of ~ 100 m/s for cool models, to ~ 600 m/s for hot or low logg 
models. Our M-star structure models are calculated with the 3D radiative-hydrodynamics (RHD) 
code CO^BOLD. The spectral synthesis on these models is performed with the line synthesis code 
LINF0R3D. 
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INTRODUCTION 

The measurement of line broadening in cool stars is in general a difficult task. For 
example, the investigation of the rotation- activity connection among field M-dwarfs 
requires the measurement of rotational line broadening with an accuracy of 1 km/s 
The spectral lines have to be very narrow and well isolated to detect slow rotation. 
But in cool stars most of the individual atomic lines become very weak at these low 
temperatures and are dominated by pressure broadening. Also the measurement of the 
magnetic field strength is dependent on the line width, and detection of Zeeman splitting 
becomes more difficult at low temperatures due to the aforementioned reasons. Since it is 
possible to use the narrow and well isolated FeH molecule lines in cool stars to determine 
radial velocities or magnetic field strength [2], it would be very helpful to characterize 
the impact of macroscopic velocity fields - primarily driven by convection - on the 
line shapes from hydrodynamical model atmospheres. We calculate 3D-C0^B0LD stellar 
atmosphere models jst] which serve as input to the line formation program LINF0R3D 
[based on 3]. Velocity fields and thermal inhomogeneity are naturally represented in 
the hydrodynamical models. The influence on the modeled spectral lines can then be 
investigated and translated into effective micro- and macro-turbulent velocities used 



in classical analyses based on ID hydrostatic model atmospheres. The comparison 
with ID-models provides some insight when it is necessary to apply 3D-models in the 
spectral analysis of cool stars. 

In the first part of the paper we describe our hydrodynamical model atmosphere code 
and give an overview of the models, in the second part we investigate the velocity fields 
in the models and their dependence on log^ and Teff. We will express the influence of 
velocity fields on line shapes in terms of micro- and macro-turbulent velocities. 



METHODS AND MODELS 

CO^BOLD is the abbreviation for "Conservative COde for the Computation of Com- 
pressible convection in a BOx of L Dimensions with L=2,3" [i5J. It can be used to 
model solar and stellar surface convection. In solar-like stars, a CO^BOLD model repre- 
sents the 3D flow geometry and its temporal evolution in a small (relative to the star's 
radius) Cartesian domain at the stellar surface ("box in a star" set-up). The spatial size of 
the domain is chosen to be sufficient to include the dominant convective scales, i.e. the 
computational box is large enough to include a number of granular cells at any instant in 
time. A CO^BOLD model provides a statistical realization of the convective flow. In this 
investigation we usually average over five flow fields taken at different instances in time 
("snapshots") to improve the statistics. 

CO^BOLD solves the coupled non-linear time-dependent equations of compressible 
hydrodynamics coupled to the radiative transfer equation in an external gravitational 
field in 3 spatial dimensions. As set of independent quantities are chosen the mass 
density p, the three spatial velocities v^, Vy, and v^, and the internal energy With these 
quantities, the 3D hydrodynamics equations, including source terms due to gravity, are 
the mass conservation equation 

dp ^ dpv^ ^ dpvy ^ dpv, ^ ^ 
dt dpx dpy dpz 

the momentum equation 

xVx+i'N ^ / pVxVy \ ^ I pVxV^ 

^1 \ ^3'x \ ''"'^^ +^ pVyVy + P +^ pVyV, 1 = 1 pgy I , (2) 

and the energy equation which includes the radiative heating term Q^ad 

H " ^ 'd'z ^ ^ 

£ik denotes the sum of internal and kinetic energy. The gas pressure P is related to den- 
sity p and internal energy via a (tabulated) equation of state P — P{p,£i). For the 



local models used here the gravity field is given by the constant vector ^ = I 




CO BOLD uses the convention that the vertical axis points upwards. The radiative heat- 
ing term Q^ad is obtained from the solution of the non-local frequency-dependent radia- 
tive transfer equation. The frequency dependence of the radiation field is captured by 
considering a small number of representative wavelength bands ["opacity binning", see 
[1,0]. The resulting 3D radiative-hydrodynamic (RHD) models treat convection from 
basic physical principles and avoid approximations like mixing-length theory. In the 
following, the three dimensional data cubes of the CO^BOLD-models will be called "3D- 
models", and the spectral lines computed from three dimensional atmosphere models 
"3D-lines". 

We will characterize the velocity fields in the 3D RHD models and analyze their 
influence of FeH lines. For the ID spectral synthesis the 3D velocity fields will be 
expressed in terms of micro- and macro-turbulent velocity. For this investigation, we 
choose a set of CO^BOLD-models with Teff = 2500 K - 4000 K and log^ = 3.0 - 5.0 
[cgs]. TablefUgives the model parameters. The opacities used in all models are obtained 

TABLE 1. Overview of different model quantities. We simulated main sequence stars and 
varied logg slightly for models with changing T^g. The models with almost same T^ff were 
started at the same Tgff of 3300 K, but they converge at slightly lower or higher Teff. 



Model code 



Dim. 



Size(x,y,z) [km] Opacities Jeff [K] logg [cgs] 



d3t33j 


|30mm00wl 


3 


85000 X 85000 x 58350 


PHOENIX 


3240 


3.0 


d3t33j 


j35mm00wl 


3 


28000 X 28000 X 11500 


PHOENIX 


3270 


3.5 


d3t33j 


T40mm00wl 


3 


7750 X 7750 X 1850 


PHOENIX 


3315 


4.0 


d3t33j 


^SOmmOOwl 


3 


600 X 600 X 260 


PHOENIX 


3275 


5.0 


d3t40j 


J45mm00n01 


3 


4700 X 4700 X 1150 


PHOENIX 


4000 


4.5 


d3t38j 


|49mm00wl 


3 


1900 X 1900x420 


PHOENIX 


3820 


4.9 


d3t35j 


jSOmmOOwl 


3 


1070 X 1070x290 


PHOENIX 


3380 


5.0 


d3t28j 


jSOmmOOwl 


3 


370 X 370 X 270 


PHOENIX 


2800 


5.0 


d3t25j 


?50mm00wl 


3 


240 X 240 X 170 


PHOENIX 


2575 


5.0 



from the PHOENIX stellar atmosphere package 0] assuming an abundance mixture 
according Asplund et al. The opacity tables were computed after Ferguson et al. 
H and Freytag et al. 

In order to compare ID model atmospheres with the hydrodynamical 3D-models, we 
average the 3D-models over surfaces of equal optical depth. We obtain so-called <3D>- 
models which have the same mean thermal profile as in the 3D-models but without 
cold and hot regions which stem from the convective granulation pattern. In the <3D>- 
models the hydrodynamic velocity field is not considered. The line-broadening is treated 
in the classical way by adding isotropic Gaussian micro- and macro-turbulence. We will 
call the spectral lines of these models < 3D > -lines. 



HYDRODYNAMICAL VELOCITY FIELDS 

Spectral lines are broadened by velocity fields where the wavelength of absorption 
or emission of a particle is shifted due to its motion in the gas. Here we are mostly 
concerned with the macroscopic, hydrodynamic motions but have in mind that the 
thermal motions are also constituting an important contribution. If we envision each 
voxel in the RHD model cube to form its own spectral line, the whole line consists of 



a (weighted) sum of single lines. The velocity distribution might be represented by a 
histogram of the velocities of the voxels. We try to describe the velocity fields in that 
sense instead of using the rms- velocities. Since a velocity vector is assigned to each 
voxel, we apply a binning method and count all vertical velocities to plot them in a 
histogram with a bin size of 25 m/s. We took the standard deviation C7 of a fitted Gaussian 
normal distribution as a measure for the velocity dispersion in the models (see FigH]). At 
this stage of the investigation time we only concentrate on the vertical direction which 
is appropriate if we assume that the major part of the line broadening stems from the 
vertical motions. The velocities of the horizontal directions are of roughly the same 
order. As a measure for the statistical variability in the 3D RHD model, we average 
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FIGURE 1. Histogram of the velocity distribution in vertical direction. The normalized number of 
points is plotted against the vertical velocity in m/s (solid line). The Gaussian (dashed line) fits the 
velocity distribution and determined an FWHM value (dashed-dotted line), which is related to a with 
FWHM — 2v21n2 • a. The underlying model is located at Teff = 2800 K and \ogg = 5 [cgs]. 

over five different temporal windows and also compute the standard deviation. The 
determined velocities are plotted against Tgff and logg to investigate their dependence 
(Fig. (21). The velocities increase strongly with decreasing log^ and with increasing 
Teff. They are as low as a few hundred m/s in the coolest high-gravity models. We 
can describe the dependence of the velocity dispersions on Tgff with a second order 
polynomial, and in \ogg with a linear function (The fits for the velocity dispersion are 
overplotted in Fig 111)- 

In order to investigate in which region of the model atmospheres the velocity fields are 
generated and what the physical meaning of the o is, we bin the 3D model velocities 
at constant geometrical height and fit a Gaussian to determine the velocity dispersion 
for each layer (Fig. [3]). Typically, we find maxima of the velocity dispersion around 
logT = — 1. The maximum velocities in the convection zone of each model (Teff 
and log^) are plotted in Fig. |2]too. In order to measure the velocities in the region 
where the lines originate, we compute the mean of the velocities weighted with the 



contribution function of the line depression [|lll] o^eighted 



showed in the 



bottom panel of Fig. [3l These velocities are also plotted in Fig. [2l In most cases they lie 
between the velocity dispersion and the maximum velocity in the convection zone. Some 
models show a strong increase of the vertical velocity in higher layers (Fig. [3]). This is 
related to the presence of waves which are excited by the stochastic fluid motions ^ 
and references therein]. However, it will not affect the spectral lines, because they are 
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FIGURE 2. Plotted are velocity dispersions with error bars, maximal velocities in the convection zone, 
and weighted velocity dispersions (see text) for models with different T^fs (left) and different logg (right). 
The dashed line on the left plot shows a second order polynomial fit for the velocity dispersion and on the 
right plot a linear function is sufficient to fit the log^ data. 
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FIGURE 3. Upper panel: a velocities. In each upper panel the a -velocities (solid) and the mean a- 
velocity (dotted) are plotted against the optical depth in logarithmic scale. Bottom panel: In each bottom 
panel are the contribution functions (CF) of a FeH-line at a wavelength of 9954.0 AThe width (dashed 
triple dotted) and depression (dashed) of the line as a function of optical depth on a logarithmic scale. 



generated in the region between an optical depth of log T = 1 .0 and log T = — 4.0. 
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FIGURE 4. Micro-turbulent velocities as a function of log gf for different Jeff (left) and different logg 
(right). 
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FIGURE 5. FeH lines for models with Tiff ^ 2800 K, logg = 5.0 [cgs] (left) and T^g = 3820 K, 
logg — 4.9 [cgs] (right). The upper panels show the 3D-line (dots) and the <3D> comohed-l^^^ (solid 
line) which was broadened by a Gaussian profile. For comparison we plotted a <3D>-line which was not 
broadened by any velocities (dashed line). In the lower panels are the 3D-<3D>„,„v„/,,grf residuals plotted. 
One can see the asymmetry which stems from the line shifts due to convective motions. 




FIGURE 6. Micro- and macro-turbulent velocities as a function of T^-g (left) and logg (right). The data 
points are fitted by second order polynomials. 



MICRO- AND MACRO-TURBULENT VELOCITIES 



In the following we express the 3D hydrodynamical velocities in terms of classical 
micro- and macro-turbulent velocities [see e.g. [l^, Ii3]. We try to see how accurately 
the ID broadening resembles the 3D broadening in the regime of cool stars which, as 
we saw, exhibit velocity fields of relatively small amplitudes. AID treatment would save 
a lot of CPU-time, and little differences between ID and 3D broadening could favor the 
usage of fast ID atmosphere codes to simulate M-stars for comparisons with observa- 
tions, e.g. to determine rotational- or Zeeman-broadening. 

Micro-turbulent velocities: In order to determine the micro-turbulent velocities, we use 
the curve of growth (CoG) method [13]. We increase the line strength of an FeH and 
an Fel line by (artificially) increasing the \o%gf value of the line which results in an 
increasingly saturated lines. In saturated lines, the effective absorbing cross-section is 
enhanced due to the influence of the micro-turbulent velocities. From 3D spectral syn- 
thesis we obtain a 3D CoG, which can we compare with the <3D> CoG for different 
micro-turbulent velocities which are included in the ID spectral synthesis. We use a 
grid of micro-turbulent velocities that range from km/s to 1 km/s in 0.125 km/s steps 
and obtain <3D> CoGs for nine different micro-turbulent velocities. We compare the 
<3D> CoGs with the 3D CoG and select the velocity of the <3D> CoG that fits the 3D 
CoG best. The determined velocities are plotted in Fig|6l Since strong (saturated) lines 
tend to be formed in the upper atmosphere, and weak (unsaturated) lines in deeper lay- 
ers of the atmosphere, \oggf and the height of formation are related. If we fit a micro- 
turbulent velocity to each \oggf point of the 3D CoG, we obtain a height-dependent 
velocity structure (see Fig|4]). 

Macro-turbulent velocities: For modeling macro-turbulent broadening we used the 
radial-tangential profile from Gray [14], as well as a Gaussian profile in the ID spectral 
synthesis. We compared <3D> absorption lines which were broadened with both pro- 
files, and the results differ only very little. Hence it is appropriate to use the Gaussian 
profile here. It is remarkable that the broadening due to small 3D RHD velocity fields 
in M-stars can be approximated with a simple Gaussian velocity distribution. Two ex- 
amples of Gaussian broadened < 3D > FeH lines are shown in Fig. [51 To determine the 
macro-turbulent velocities of the 3D RHD models, we first compute the < 3D > lines 
with a given micro-turbulent velocity and after this we broaden them with a Gaussian 
with a given macro-turbulent velocity until they match the 3D profiles. The results are 
plotted in Fig|6l The broadened < 3D > FeH lines fit the 3D FeH lines very well. The 
difference of the < 3D > and 3D centroid (C = ^^) range in the order of a few m/s for 
small velocity fields to 30 — 40 m/s for strong velocity fields in hot M star models, or 
models with low logg. The differences stem from the asymmetry of the 3D line profiles, 
which are shifted due to convective motions. However, the errors in the normalized flux 
profiles are lower than 1% (see Fig. [5]). We fit the micro- and macro-turbulent veloc- 
ities as a function of T^ff and log^, each pair can be fitted with a polynomial second 
order. The fits are overplotted in Fig [61 Both micro- and macro-turbulent velocities show 
similar dependence on Teff and logg as the RHD velocity fields considered before (i.e. 
they increase with increasing Tgff and decreasing logg). Especially the macro-turbulent 
velocities for changing T^ff agree well with the determined RHD velocities. 



SUMMARY 



We investigated the velocity fields in a set of M-star 3D hydrodynamical models. They 
range in Teff and log^ from Teff = 2500 K - 4000 K and logg = 3.0 - 5.0 [cgs]. 
We applied a binning method to characterize the velocity structure and to determine 
mean velocities. These velocities range from ~ 100 m/s for the cool models up to 
~ 800 — 1000 m/s for hot models or models with small logg values. 
In order to compare the 3D RHD velocity fields with velocities needed in ID spectral 
synthesis, we expressed them in terms of classical micro- and macro-turbulent velocities. 
For this purpose we computed the 3D and <3D> Curve of Growth and determined 
the micro-turbulent velocity from them. They range in the order of ~ 100 m/s for cool 
models to ~ 600 m/s for hot models or models with small logg. The macro-turbulent 
velocities were determined through convolution with a Gaussian profile which, what 
turns out, is sufficient in the regime of cool stars. The obtained velocities are of the 
order of the 3D RHD velocity fields and show similar dependence from Tgff and logg. 
We saw that ID spectral synthesis in cool stars with micro- and macro-turbulent velocity 
broadening is able to reproduce the 3D spectral line synthesis. 
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